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Abstract. For an ideal I with a positive dimensional real variety, based on 
moment relaxations, we study how to compute a Pommaret basis which is 
simultaneously a Groebner basis of an ideal J generated by the kernel of a 
truncated moment matrix and nesting between I and its real radical ideal. We 
provide a certificate consisting of a condition on coranks of moment matrices 
for terminating the algorithm. For a generic delta-regular coordinate system, 
we prove that the condition is satisfiable in a large enough order of moment 
relaxations. 



1. INTRODUCTION 

Finding real solutions of a polynomial system is a classical mathematical prob- 
lem with wide applications. Let / = {hi, . . . , hm) C W[x] :— R[xi, . . . , a;„] be an 
ideal generated by polynomials hi, ... , h„i € R[x]. Its complex and real algebraic 
varieties are defined as 

Veil) {x e C" I fix) = V/ e /}, Vm{I) := Vc{I) D W\ 

The vanishing ideal of a set y C C" is an ideal 

I{V) := {/ e C[x] I f{v) = 0,VveV}. 

The radical (also called complex radical) of / is 

V7 {/ e C[x] \ f'' e / for some fc G N} , 

while the real radical of / is defined as 

^ := |/ e R[x] +J2qf el for some /c G N, gi, . . . , G M[a;]| . 

Clearly, they satisfy the inclusion / C ^/I C v^. An ideal / is called radical (resp. 
real radical) if / = ^/l (resp. / = V^). According to the Real Nullstellensatz 
0, the vanishing ideal /(Vr(/)) of the zero set Vr(/) is a real radical ideal and 
I{Vu{I)) = ^. 

There exists much work on computing a complex radical ideal vT, like [3 El [TTJ 
[T^. The algorithms range from numerical ones (e.g., [131 [THl [H]) to symbolic 

ones (e.g., [6l|3T]). For the general case of / being positive dimensional, a commonly 

used technique is to reduce the problem to the zero-dimensional case, like in Gianni 

et al. [n] and Krick and Logar [l4] . 

The problem of computing the real radical ideal is typically much more 

difficult than computing -\/T. Becker and Neuhaus [1] proposed a symbolic algorithm 
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based on the primary decomposition to compute (also see [351 [Ml [311 |3D] ) . Some 
interesting algorithms based on critical point methods were proposed in [I] [H ISl HH] 
to compute a point on each semi-algebraically connected component of real algebraic 
varieties. 

A new approach based on moment relaxations has been proposed by Lasserre 
et al. [m [TH m [H] for computing v^, provided it is a zero-dimensional variety. 
Hereby we briefly describe this interesting approach. 

For a sequence y = {ya)aeji^ G , its moment matrix 

is a real symmetric matrix whose rows and columns are indexed by the set T" := 
{x" I a € N"} of monomials. Given a polynomial h G R[x], set vec(/i) := (/lajaeN- 
and define the sequence hy := M{y)vec{h) G R^". We say that a polynomial p hes 
in the kernel of M(y) when M{y)p := M{y)vec{p) = 0. Given a truncated moment 
sequence y = {ya)aew^^ & M^^t^ it defines a truncated moment matrix 

Mt{y) ■■= {ya+l3)a..l3mi' 

indexed by the set := {x" | a G NJ* with |a| := Sf^itti < t}. 

We work with the space M\x\t of polynomials of the degree smaller than or equal 
to t. For a polynomial p £ M[a;]t, if Mt{y)vec{p) ~ 0, we say p lies in the kernel of 
Mi(y), i.e., 

(1) kerMt(y) := {p e R[x]t \ Mt{y)vec{p) = 0}. 
Let / — {hi, . . . , h„i) C M.[x] be an ideal and set 

(2) dj rdeg(/ij)/2], d := max dj. 

l<.j<m 

For t > define the set 

(3) ICt {y e | yo = 1, Mt(y) h 0, Mt-d, {h,y) = 0, j = 1, . . . , m}. 
An element y ICt is generic if Mt{y) has maximum rank over JCt- We denote 

(4) /Cj*^" {y (z JCt \ rankAft(?/) is maximum over JCt}. 

When the real algebraic variety Vffi(/) is finite, Lasserre et al.[T7] used the flat 
extension (a rank condition of moment matrices in |^) as a certificate to check 
whether polynomials in ker Ms{y) (1 < s < for a generic element y £ JCt generates 
the real radical ideal /(Vk(/)). When Vk(/) is positive dimensional, this certificate 
does not work. The example given by Fialkow in [TDj Example 3.2] can be used 
to explain the difficulty. Unlike the zero-dimensional case, although the kernel of 
the moment matrix of the third order consists of only a polynomial z ^ which 
is already a Grobner basis of the real radical ideal / = /(Vr(J)) = {z — x^), we can 
not extend the truncated moment sequence y G JC3 to the next order, i.e., y has no 
representing measure. 

The motivation of this paper is to provide a certificate for checking (kei Mt{y)) = 
J(Vr(/)) when Vr{I) is positive dimensional. Unfortunately, we still can not solve 
this open problem [22] §2.4.3] completely. However, we provide a certificate PU)) 
based on the geometric involutivity theory [321 1301 133 for checking whether we 
have obtained a weak Pommaret basis (also a Grobner basis) of an ideal J = 
(kerMt_2(y)) satisfying / C J C /(Vm(/)) under graded reverse lexicographic or- 
der. A (weak) Pommaret basis is a special form of the familiar Grobner basis 
which allows for directly reading off the depth, the projective dimension and the 



A CERTIFICATE FOR SEMIDEFINITE RELAXATIONS IN COMPUTING POSITIVE DIMENSIONAL REAL VARIETIES 

Castelnuovo-Mumford regularity of a module. When the real algebraic variety Vs.{I) 
is positive dimensional, for examples in Section |4l we succeed in showing that the 
computed Pommaret basis is an involutive basis of the real radical ideal /(Vk(J)). 
In general, it is still not possible to prove that the kernel of the moment matrix 
satisfying the certificate ((T0| generates a real radical ideal. 

The paper is organized as follows. In Section [21 we review some preliminary 
backgrounds like elementary algebraic geometry, moment matrices, involutive di- 
visions and involutive bases. In Section [31 we present an algorithm based on the 
semidefinite programming and moment relaxations in computing a Pommaret basis 
of an ideal J satisfying / C J C Vm{I) = Vc{J) n M", and propose a cer- 

tificate for terminating the algorithm and prove it works for a positive dimensional 
Mr (7) under a 5-regular coordinate system. In Section [H we present computational 
results for a set of examples in [371 1301 IS21 [31] • Some open questions and ongoing 
work are given in Section [5l 

2. PRELIMINARY 

We introduce some notations and preliminaries about polynomials, matrices, 
semidefinite programs and the involution. Given IK = K or C, the ring of multivari- 
ate polynomials in n variables over the field K is denoted by K[a;] :— K[a;i, . . . , a;„]. 
For an integer t > 0, K[x]t denotes the set of polynomials of degree at most t. N de- 
notes the set of nonnegative integers and we set N" := {a G N" | \a\ := E^^j^a^ < t} 
for t e N. For a £ N", x" denotes the monomial x"^ ■ ■ ■ x"" whose total de- 
gree is \a\ := T^f^iUi. All monomials are included in T" := {x" | a G N"} and 
T" :— {x" I a S N"} consists of monomials with degrees bounded by t G N. Con- 
sider a polynomial p € K[a;], p — T,a<£N"PaX°' , where there are only finitely many 
nonzero pa € K, its leading term It^(p) is the maximum term x" with respect to a 
monomial order ^ for which p^ ^ 0. We denote by (lt^(/)) the ideal generated by 
leading terms of polynomials in /. The symbol [x]t denotes the sequence consisting 
of all monomials of degrees less than or equal to t: 

[x]f :~ [1, Xi , ■ ■ • , X„, X]^,XiX2,''' ,X]^,X-^ X2, ■ ■ ■ , x„] . 

2.1. Properties of moment matrix. The kernel of a moment matrix is particu- 
larly useful as it has the following properties, see [8l [TZl [20l [2T1 [24] . 

Lemma 2.1. [17l Proposition 3.6] Let]ierM{y) := {p e M[x] | M{y)vec{p) = 0} be 
the kernel of a moment matrix M{y). Then ker M{y) is an ideal in R[x] . Moreover, 
if Mijj) >: 0, then kerM(j/) is a real radical ideal. 

The kernel of the truncated moment matrix Mtiy) is not an ideal, but under 
certain conditions, it has the following properties. 

Proposition 2.2. [TT! Lemma 3.5, 3.9] Let y G K^?* and its truncated moment 
matrix Mt{y) is positive semidefinite. 

(i) Iff,ge M[x] with deg(/g) <t-l, then f G kei Mtiy) =^fge kerMt(y). 

(ii) For a polynomial p G K[x], if p^^ + cr G Vev Mt{y) for some fc G N and 
ctgX^RH^. thenpGkeiMt{y). 

(iii) We have kerM((y) n K[x]s = ker Ms{y) forl<s<t. 

Generic elements of /Ct have useful properties. The following results are cited 
from [m Lemma 3.1] and [27l Lemma 7.28, 7.39]. 
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Proposition 2.3. Assume y € ICf is generic. 

(i) For all 1 < s < t, we have kei Ms{y) C \/l and 'ker Ms {y) C kerMs(z) for 
all z G /Cf . 

(ii) Ift<t' and y' e /C^™, then ker Mt(y) C kerMt-(y')- 

(iii) For every finite basis {gi, . . . ,gk} of the real radical ideal \fl , there exists 
to G N such that gi, . . . , gk € ker Mt(z) for all z d JCt and t > to- 

(iv) It holds that (kei AIt{y)) = \fl if t is sufficiently large. 

In the following, we review some properties of moment matrices in the occurrence 
of inequality constraints. Consider the semialgebraic set 

(5) A := {x e R" I f,{x) > 0, . . . , fsix) > 0}, 

where /i, . . . , G R[a;]. The ^-variety Va{I) denotes the intersection 

VAii)^VM{i)nA. 

For every e {0, 1}'*, we denote the product f :— /j;'V2^ • ■ ' fs°- 
Definition 2.4. |23) The A-radical of an ideal I is defined as 



p e R[x\ 



p^'' + a^f e / for some keN,a^eY^ R[xY 



!^e{o,i}'' 

The ideal I is called A-radical if I = \fl . 

Theorem 2.5. [551 Semialgebraic Nullstellensatz] Let I he an ideal in M.[x] and A 
be defined by (0)- Then \fl is an A-radical ideal and \fl = /(Vr(/) H A). 

To compute the yl-radical ideal \fl^ we consider the set 

(6) Kt,A ■■= ICtn{ye R^- : Mt„rf^.„ {fy) h 0, Vi. G {0, 1}^} , 

where df I' — [deg(/'')/2] . Clearly, the set /Ct^^ is a restriction of /Ct. The definition 
of the set JCt^A is motivated by the polynomials in and the Semialgebraic Null- 
stellensatz. The generic elements of JCt^A are similarly defined to be the elements 
of the set 

/Cj^ := {y e ICt^A ■ rankMt(?;) is maximum over ICt.A}- 

Lemma 2.6. Let {gi, . . . ,gk} be a set of generators for the ideal yfl. Then there 
exists to ^ N such that gi, . . . ,gk € ker Mt{y) for all y G ICt.A and t >to. 

The following proof mimics the proof of Claim 4.7 in [T7] . 



Proof. For each £ — 1, . . . , fc, by Theorem [231 there exists to; e N and polynomials 
a^, e SR[a;]^ and uj € R[x] for 1 < j < to such that 

ni 

(7) gf"^'+ J2 ^^r = Y.^,h,. 

i/e{o,i}= J=i 

For t > to, where 

to = l + ina.x{d,deg{gf"''),deg{a^f''),deg{ujhj)), 
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since deg{ujhj) <t — l and hj e ker Mt{y), by Proposition [221 (i), we have ujhj E 
keiMtiy), i.e., gf'^' + E.e{o,iy '^•^f e keiMtiy). Set a, = e SR[x]2, 

then we have 

vec(gr')^^t(y) vec(gr') + ^ vec{a,^,f Mt^df. [Fy) vec(a,,,) = 0. 

Since Mt{y) >z. and Mt-dfuif^y) ^ 0, every summand in the above expression 
must be zero, and thus 5™' € ker Aft(j/). If m; is even, g™' e kerMt(?/) imphes 
^m,/2 ^ kerM((T/). If m; is odd, since deg((7™'^^) < t — 1, we have 

5,"' e kerMt(y) ^ .9;"'+^ e ker Mt{y) ^ g^^'+^'Z^ g kerMt(?;). 
Repeat this process, we can show that gi e ker Aft(y). □ 

Theorem 2.7. There exists io e N such that (ker Aft (y)) = \/T /or y e /Cf^' 
and t > tQ. 

Proof. Let {g'l, . . . , (7^} be a set of generators for the ideal According to Lemma 
12. 6[ we can choose to £ N such that gi, . . . ,gk € ker ^^((y) for all y £ JCt,A and 
t > to. Let y G /Cf™ and choose an arbitrary point v G Then [v]2t € /Cj,^ 

and z = (y + [w]2t)/2 £ /C*,^. Clearly, it holds that 

ker Aft((y + [v]2t)/2) = ker Aft(y) n ker Aft(H2t). 

The rank of Mt{y) being maximum implies ker Aft((y + [w]2t)/2) — ker Mt{y) and 
ker Mt{y) C ker Aft([?;]2t). For every p G ker Aft (y), we must have p G ker Aft ( [w] 2t ) 
and p{v) = 0. This means that p vanishes on the set Va{I)- By Theorem 12.51 we 
get p G \fl and thus the inclusion (kerAft(y)) Q \fl = I{Va{I)) holds. Since 
5i,...,yfc G kerAft(y), we get (kerAft(y)) D ^ = I{Va{I)) and the proof is 
completed. □ 

2.2. Involutive Divisions and Involutive Bases. When the real algebraic va- 
riety T/R(f ) is finite, Lasserre et al. [16l [17] proposed new approaches based on 
moment relaxations for computing Grobner bases or border bases of the real radi- 
cal ideal For the positive dimensional real variety 14 (f), we can also compute 
its Grobner bases. Stimulated by the work in [18] and [26j [29j [30] , we propose a new 
approach based on the completion to involution to compute a Pommaret basis of 
an ideal nested between f and v^. A Pommaret basis is simultaneously a Grobner 
basis, but contains extra information such as the Castelnuovo-Mumford regularity. 
Moreover, we provide a new stopping criterion for the algorithm which is based on 
the classical Cartan's test for involution from the theory of exterior differential sys- 
tems. We now introduce some basic concepts from the classical theory of involutive 
systems for polynomial systems. For background, see |32i 133) . 

Definition 2.8. Let v — [z^i, . . . , G N" he the multi index of a monomial x" . If 
k is the smallest value such that Vk ^ 0, then the class of v or is k, written by 
cls(z^) — k or cls(x'^) — k. The class of a polynomial f which is denoted by cls(/) 
is k, if the class of its leading term cls(lt^(/)) = k. 

We say that a term order respects classes, if for monomials x'^ and x'' of the 
same total degree, cls(/i) < cls(z^) implies x^ -< x'^ . An important example of a 
class respecting ordering is the graded reverse lexicographic order ^tdeg. 
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Definition 2.9. With an ordering on the variables xi -< ■ ■ ■ -< Xn, the graded 
reverse lexicographic order ^tdcg is defined by x" ^tdcg x^ , if \a\ < |/3|, or \a\ = 
and the first non-vanishing entry of the multi index a — /3 is positive. 

Throughout the paper, we use ^tdog in assigning orders of monomials, and sorting 
rows and columns of a moment matrix Mt{y). Let (N",+) be an Abelian monoid 
with the addition defined componentwise. For any multi index S N" , we introduce 
its cone C(y) = v -\-W^, i.e., the set of all multi indices that can be reached from v 
by adding another multi index. 

Definition 2.10. [331 Definition 3.1.1] An involutive division L is defined on the 
monoid (N",+), if for any finite subset B C N", a set N]^t}{v) C {l,...,n} of 
multiplicative indices, and consequently a submonoid L{u,B) = {/i e N" | Vj ^ 
Nl,b{v) ■ fJ-j = 0}, is associated to every multi index v € B such that the following 
two conditions on the involutive cones Cl.b(i^) = h' + L[v,B) C N" are satisfied. 

(i) // there exist two elements ^ B with CL.BifJ-) (~^Cl,b{i^) 9, either 
CL.Biti) Q Cl,b{^') orCL.Bi^) C CL.Bilj) holds. 

(ii) If B' C B, then Nl,b{v) C Nl.,b'{v) for all v G B' . 

An arbitrary multi index /i G N" is involutively divisible by v € B, written v \l,B fJ-, 
if fJ- € Cl,b{v). In this case v is called an involutive divisor of fi. 

Definition 2.11. |331 Example 3.1.7] The Pommaret division written by P is de- 
fined by a simple rule: if els (v) = k, then we set Nl,^{i') — {1, . . . , k}. 

Remark 2.12. The Pommaret division is a globally defined division as the assign- 
ment of the multiplicative indices to a multi index ly G B is independent of the set 
B. The Pommaret division is an involutive division by [33i Lemma 3.1.8]. 

Definition 2.13. [33, Definition 3.1.9] The involutive span of a finite set B C N" 
is 



The set B is called weakly involutive for the division L or a weak involutive basis 
of the monoid ideal {B) , if (B) l={B) . The set B is a strong involutive basis or for 
short an involutive basis, if the union is disjoint, i.e., the intersections of the 
involutive cones are empty. 

For a polynomial / G K[x] and a term order ^, we select its leading term 
lt^(/) = x^ with the leading exponent le^(/) = /i. 

Definition 2.14. [33, Definition 3.4.1] Let I C K[x] be an ideal. A finite set H C I 
is a weak involutive basis of I for an involutive division L on N", if \e^{'H) is a 
weak involutive basis of the monoid ideal le^(/). The set % is a strong involutive 
basis of I, if\e^{'H) is a strong involutive basis o/le^(/) and two distinct elements 
of % never possess the same leading exponents. 

Remark 2.15. Definition 12.131 and Definition 12.141 implv immediately that any 
weak involutive basis is a Grobner basis. 

Not every ideal in K[a;] possesses a finite Pommaret basis (see [33]). 

Definition 2.16. [33t Definition 4.3.1] A coordinate system is called 5 -regular for 
the ideal I C 'K.[x\ and the term order -<, if I possesses a finite Pommaret basis for 
the term order ^. 
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Theorem 2.17. [33) Theorem 4.3.15] Every polynomial ideal I C K[a::] possesses a 
finite Pommaret basis for a term order -< in suitably chosen coordinate systems. 

Definition 2.18. [531 Definition 3.4.2] Let T C ]K[a;]\{0} be a finite set of polyno- 
mials and L be an involutive division on N". We assign to each element f (z J- a 
set of multiplicative variables 

The involutive span of J- is then the set 

Theorem 2.19. [33l Theorem 3.4.4] Let L C K[a;] be a nonzero ideal, U C /\{0} 
a finite set and L an involutive division on N". Then the following two statements 
are equivalent. 

(i) The set % G I is a weak involutive basis of I with respect to L and -< . 

(ii) Every polynomial f £ I can be written in the form 

(9) f^Y.p^-^ 

hen 

with coefficients Pt G IK[Xi^-^.^(/i)] satisfying lt^(P/j • h) ^ \i^{f) for all 
polynomials h Cz H such that Ph =/= 0. 

H is a strong involutive basis, if and only if the representation (0) is unique. 

Corollary 2.20. [33l Corollary 3.4.5] Let H be a weak involutive basis of the 
ideal I C K[x]. Then {T-Oj^^^ — L. If % is even a strong involutive basis of 
I, then I considered as a K-linear space possesses a direct sum decomposition 

Proposition 2.21. [33l Proposition 3.4.7] Let I C K[x] be an ideal and H C I be 
a weak involutive basis of I for the involutive division L. Then there exists a subset 
% (1% which is a strong involutive basis of L 

Definition 2.22. // we regard K.[x\ as a linear space, then the ideal I and the 
truncated ideal It = I C] K[x]( are both subspaces in K[a;]. We say that the set 
G — {gi, . . . ,gs\ is a reduced basis of It, if it is a linear independent basis of It and 
all polynomials in G have different leading monomials with respect to a given term 
order. 

3. COMPUTING A POMMARET BASIS 

In this section, we present an algorithm as well as a certificate for computing 
a Pommaret basis for an ideal J, s.t. / C J C when Vr(/) is positive 

dimensional. The certificate given in (jlOp generalizes the flat extension conditions 
in [17] for the zero-dimensional real variety to the positive dimensional case. 
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3.1. The certificate. Let / = {hi, . . . , h,n) C M.[x] be an ideal and d := inaxi<j<m dj, 
dj \deg{hj)/2]. For each t > d, recall the notions 

ICt:^{ye Kf*?' I yo = 1, Mt(y) h 0, Mt_d, {h.y) ^ 0,j = 1, . . . ,m}, 

and 

j^gen ^ jQ^ I rankMt(t/) is maximum over lCt\- 

For a moment matrix Mt{y) of order t, the truncated moment matrix Mt-i{y) 
for I <t IS the order t — £ principal submatrix of Mt{y) indexed by a, /3 G N"_^. 

Let ckj denote the number of class j polynomials of degree i — 2 in the reduced 
basis of ker Aft_2(?/). We have the following theorem: 

Theorem 3.1. Let I — {hi, . . . , hm) C W[x] be an ideal. If there exists an integer 
t > 2d satisfying 

n 

(10) f""^ keiMt-2{y) = corankMt_i(2/) - corankMt_2(y), 

for an element y € /C^*^". Then a reduced basis of the null space of Mt-2{y) is a 
weak Pommaret basis for J = (ker Af(_2(?/)) under the monomial ordering ^tdcg 
and 

(11) ^ C (J) C I{Vm{I)), Vm{I) = VciJ) n M". 

The proof of Theorem 13.11 follows from Proposition 12.31 and Theorem 13.91 whose 
proofs are given in Section [3731 

Remark 3.2. Although the reduced bases of ker Mt-2{y) are not unique, they have 
the same set of leading terms since they can be represented linearly by each other. 
Therefore, each reduced basis of kerMt_2(y) has the same value of J2]=ij'^j- 

In our algorithm, we need to find an element y in /Cf maximizing the rank of 
Mt{y). As pointed out in [T7], this could be done typically by solving the semidef- 
inite program 

(12) min s.t. y € /Ct 

with interior-point algorithms using self-dual embedding, see |37| I38j . 

3.2. An algorithm for computing a Pommaret basis. We list main steps of 
our algorithm based on solving (IT2|) for computing a Pommaret basis of the ideal 
J = {kci Mt-2{y)) nested between / and /(Vr(/)). 

Algorithm 3.3. Computing a Pommaret basis of an ideal J such that / C J C 
^(Vk(/)). 

Input: A set of polynomials {hi, . . . , hm} generating / and the monomial or- 
dering ^tdcg on variables xi, . . . ,Xn- 

Output: A Pommaret basis for (ker Mt_2(y)) under the monomial ordering 

^tdcg- 

Step 1: For t > 2d, compute a generic element y € /Ct by solving 

Step 2: Compute a reduced basis of kerMt_i(y). Let {gi,. . . ,gs+t} be poly- 
nomials of degree t—2 in this reduced basis. Compute the value of J2^=i j^j ■ 

Step 3: Compute corankMt_i(y) — corankMt_2(y) by calculating the num- 
ber of polynomials of degree i — 1 in the reduced basis of kerMt_i(y). 

Step 4: Test whether the condition ([T(I|) is satisfied. 
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• If yes, {gi, . . . , gs+r} is a weak Pommaret basis for (ker Mt-2{y)) and 
can be reduced further to a (strong) Pommaret basis. 

• Otherwise, let i := t + 1 and go to Step 1. 

In Section 13.31 we prove that Algorithm 13.31 is correct and terminates in a finite 
number of steps in a (5-regular coordinate system for v^. The algorithm has been 
implemented in Matlab using the GloptiPoly toolbox 12^ and we demonstrate its 
performance on a set of examples in Sectional 

Remark 3.4. In order to check the condition ([TU]). we need to compute a reduced 
basis of the null space of the truncated moment matrix Mt^i{y). These compu- 
tations have to be performed stably. For the computation of a reduced basis it 
is important to choose a proper tolerance to ensure that there is no information 
missing in kerMt_i(y). We list the tolerance used for each example in Section 4. 

3.3. Justification of the certificate. Our main goal in this section is to prove 
that Algorithm 13.31 is correct and it terminates after a finite number of steps in a 
(5-regular coordinate system for v^. 

Assumption 3.5. Let I = {hi, . . . , hm) Q M[x] be an ideal. Suppose there exists 
an integer t > 2d satisfying the condition ilO\) for y € ICf^" . Let {gi, . . . , gs+r} be 
a reduced basis of kei Mt^2{y) , where 

deg((7i) = t — 2 for 1 < i < s, and deg{gi) < t ~ 2 for s + 1 < i < s + r. 

Lemma 3.6. Under Assumvtion \3.5[ the polynomial set 

{xigi, . . .,Xj^gi, . . . ■ • .,Xj^gs,gi, . . . 

is a reduced basis o/ker Af4_i(y), where ji — c\s{gi) for i — 1, . . . , s. 

Proof. For k — I, . . . ,n, i — 1, . . . , s + r, since deg{xkgi) < ^ — 1, by Proposition [221 
(i), we have Xkgi G ker A/t_i(?/). In fact, since each polynomial in {gi, . . . ,gs+r} 
has different leading terms, according to Definition l2.101 the polynomials 

(13) xigi, . . .,Xj,gi, . . . ,xigs, . . . ,Xj^gs 

all have distinct leading terms of degree t—l. Hence they are linearly independent. 
Suppose there are aj polynomials of class j in {171, . . . , gs}, then polynomials in (jl3p 
yield j'^j hnearly independent polynomials of degree t— 1 in kerMt_i(?/). On 

the other hand, the number of linearly independent polynomials of degree t — 1 in 
a reduced basis of kei Mt-i{y) equals to corankMt_i(y) — corankMt_2(y). Hence, 
the condition (fTOll and Proposition 12.21 (iii) implies that the conclusion is true. □ 

Remark 3.7. Under Assumption 13. 5) for any polynomial / E kcrMt_i(y), we can 
express it as a linear combination; 



where Cik G M and \t^{cikXigk) ^tdeg lt^(/) for 1 < « < cls{gk) and 1 < fc < s, 
Afe € M and \t^{Xkgk) ^tdcg lt^(/) for 1 < fc < s -t- Note that every poly- 
nomial in {xigi, . . . ,Xj-^gi, . . . ,xigs, . . ■ ,Xj^gs, gi, . . ■ , gs+r} has a different leading 
term. Under the graded monomial ordering ^tdcg, there is only one Cigko 7^ with 
lt^{xiggkQ) = h^(/) if not all Cik are zeros. This property is very important and 
will be used in the proofs of theorems below. 




(14) 



k=l i=l 



k=l 
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Lemma 3.8. Under Assumption lS. 51 for all monomial and polynomials gj with 
deg(gj) < t — 2, j ^ s + 1, . . . , s + r, the polynomial x^gj can be expressed as 

s s+r 

(15) x^'gj ='^hkgk + ^ >^kgk, 

k=l k=s+l 

where hk G W[x] and Xk G M satisfying lt-j{hkgk) ^tdog lt^(2^'^gj"), fc = 1, . . . , s and 
lt^(Afc5fe) ^tdcg lt^(a;^5-,), fc = s + 1, . . . , s + r. 

Proof. If deg{x^gj) < t ~ 1, by Proposition 12.21 (i), we have x^gj G ker Aft_i(t/). 
According to Remark 13.71 we have the expression (fT5|) . Otherwise, we set x^^ = 
x^^x^"^ such that deg{x^^^ gj) = t — 1. Hence, we have 

s s+r 

x^'g,^x^'^x^-gJ^x^'^{Y.^k9k+ J2 ^'«5fe) 

fe=l k=s+l 
s s+r 

= '^x^^hkgk+ ^ XkX^^gk- 

k=l k=s+l 

We can repeat the above reduction on x^^gk for s + 1 < fc < s + r. Since deg(2:'^^ ) < 
deg(a;^), after a finite number of steps, we have the expected form (jisp . □ 

Theorem 3.9. Under A ssumvtion l 3. 51 a reduced basis {gi, . . . , gs+r} of ker Mt-2{y) 
is a weak Pommaret basis of the ideal (kerMt_2(y)). 

Proof. We show that any polynomial / G (kerMt_2(y)) can be represented as 

s s+r 

(16) f = ^hkgk+ ^ Xkgk, 

k=l k=s+l 

where G K. and hk G M[a;i, . . . , a;cis(gi.)]- Since lt^{hkgk) and lt^{gk) are all 
different for l<k<s + r,iif satisfies pB]) . then we have \t^{hkgk) ^tdeg lt^(/) 
for 1 < fc < s and lt^(Afc5fc) ^tdcg !!-;(/) for s + 1 < k < s + r. Therefore, according 
to Theorem 12.191 the polynomial set {gi, . . . ,gs^r} is a weak Pommaret basis of 
the ideal (kerMt_2(y)). 

Since {gi, . . . , gs+r} is a reduced basis of kerM4_2(j/), every polynomial / G 
(kerMt_2(j/)) can be represented as 

s+r 

where hj G R[a;], j = 1, . . . , s+r. Hence, we only need to show that each polynomial 

x'^gj for /i G N" and I < j < s + r can be written as (ITS)) . 

Set / = x^^gj. If deg(/) < t — 1, by Lemma [321 we have the expected expression 

(|16p directly. Otherwise, we prove by the induction on its leading term It^ (/) — 

to, i.e., we assume that / = x'^gj has the expected expression as long as 

lt^(/) ^tdog ^0 for G N" and l<j<s + r, we show it has the expected 

expression when lt^(/) — to- 
ll x'^ G R[xi, . . . ,XciB(gj)], nothing is to be proved. Otherwise, without loss of 

generality, let Xi^ be a non- multiplicative variable in x^ with respect to gj. Since 
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deg(gj) < t — 2, j — 1, . . . , s + r, by Proposition l2.2l fi). we have Xi^gj G ker Aft_i(y). 
By Lemma 13.61 and Remark 13.71 we have 

s cls(gfc) s+r 

k=l i=l fe=l 
s cls(gk) s+r 

(17) ^ X! X! <^ik{x^/xti)xigk + \k{x^/xi^)gk- 

k=l 2=1 fe=l 

According to Remark 13.71 there are two cases: 

(i) if all Cik = 0, there exits only one 1 < ji < s + r, such that A^^ ^ and 
lt^(Aji {x''/xi,)gjj = to; 

(ii) otherwise, there exists 1 < ji < s and 1 < ^2 < cls(gjj such that Ci^j-^ ^ 
and lt^{ci^j^[x'' /xt^)xi^gj^) = to- 

In both cases, all other terms in (I17p have leading terms of order less than tg, 
which can be expressed as (jl6|) by induction. Moreover, above two cases do not 
exist simultaneously. Therefore, we only need to check whether the polynomial 
Xj-^[x^/xi^)gj^ in case (i) or Ci^ji{x^ / Xi^^Xi^gj^ in case (ii) has the representation 

m- 

In case (i), if x^jxi^ G . . . , a;cis(gjj] then we obtain the representation 

(|16p . Otherwise, we repeat the reduction to the polynomial {x^ j Xi-^g^^. Since 
h^(Aj, (x'^/xijgjj = lt^(a;^gj) = to, we have deg(gj) < deg(5jj, i.e., 

lt^(5j) ^tdcg lt^(5ji)- 
In case (ii), if x^ jxi^ £ R[xi, . . . ,Xc\s(gj^)]i since Xi^ is a multiplicative variable 
of lt^(5jj), then (^x^ /xij^^Xi^ G M.[xi, . . . ,Xcis{gj^)]- Hence, we obtain the represen- 
tation (|16p . Otherwise, since is a non- multiplicative variable of \t-^{gj) and Xi^ 
is a multiplicative variable of It ^{gj-^), we have 

cls(gj) < cls(xij, chixi^) < cls(gjj. 
Because ^^{ci^j-^ (^x^^ / Xi^^Xi^gj^) — to, we have It^lxi^gj-^) — \t^{xi-^gj) and 

(18) cls(xi2) = cls(xi2gjj = cls(xii5j) < cls(a;jj. 

This implies that x^^ ^tdcg x,^. If lt^(gjj ^tdcg It^(gj), we have lt^(a;i2gjj ^tdeg 
lt^{xi-^gj) which leads to a contradiction. Therefore, we can deduce that 

lt^(5j) ^tdcg lt^(5ji)- 

In both cases, if the reduction does not stop, we will obtain a sequence of poly- 
nomials satisfying 

lt^(5j) ^tdcg It^(gji) ^tdeg • • • ^tdcg 1*^(3^,) ^tdog lt^(.gj,+i) ^tdog • • • ^tdcg to- 

Since the number of polynomials with strict increase leading terms bounded by 
lt^(/) = to is finite, the above procedure will stop in a finite number of steps and 
we obtain the expected form (|16p for /. □ 

Theorem 3.10. In a 6-regular coordinate system for ^fl , after a finite number of 
steps, Alaorithm \3.3\ will terminate and return an integer t > 2d which satisfies the 
condition \10\) for an element y € ICf^^^ . 
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Proof. In a (5-regular coordinate system, we have a finite Pommaret basis H = 
{hi, . . . , hs} for the real radical ideal I{Vs.{I)). According to Proposition 12.31 (iii), 
we can conclude that there exists an integer ti such that the Pommaret basis 
{hi, . . . , kg} is contained in kei Mt{y) for all y G JCt and t > ti. 

Since H is a Pommaret basis of /(Vr(/)), according to Corollary I2.20[ for t > 
ti + 2, we have the following decomposition: 

(19) I{VM{I))t-2 = ^Xi, X,i,(ft,)]t_2_deg(/,.) • hk. 

Let 

(20) T = {x"hk I e R[xi, . . . , Xeis(,,.,)], deg(x") <t~2- deg{hk), l<k<s}. 

According to Proposition!^ fi). T C ker Aft_2(y)- Therefore, by ^ and (UHl), we 
have 

/(Vi(/))t-2CkerAft_2(y). 
On the other hand, y is a generic element, by Proposition 12.31 (i), we have 

kerMt_2(y)C/(MR(/))t-2. 
Hence, we have ker Mt-2{y) = I{VM{I))t-2 and the decomposition: 

(21) kerA/4_2(y)= ^ Mxi, ■ ■ ■ ,Xcis{hk)]t-2-dcg{hk) ■ hk- 

Since "H is a Pommaret basis of /(Vh(/)), according to Definition 12.141 each 
polynomial in T has a different leading term. Therefore T is actually a reduced 
basis of ker Aff_2(2/)- By Remark 13.21 it suffices to show that the condition (fTO|) 
holds for the polynomials in T. 

Similar to the decomposition (j21[) . we can show that there exists a direct sum 
decomposition of kerMt_i(y): 

(22) kerMt_i(y) = M[xi, . . . , a::cis(/i,)]t-i-dcg(hfc) • hk- 

hken 

For a polynomial / G ker Aft_i(y) with deg(/) = i — 1, according to (P^ . we have 
the following equalities: 

s 

/ = ^ ^ c^kx^hk (note that G M[xi, . . . , Xcis(/t;,)]) 

fc=l 0<|Ai|<t-l-dog(/ifc) 
s s 

^ c^ikX^hk + ^ ^ c^kX^hk 

k=l |^i|=t-l-dcg(;ifc) fc=l 0<|/j|<t-2-dog(/ifc) 

s s 

= X! X! c^fca;cis(xf)(a;^/a;cis(2;A'))^fe + X! X! c^kX^hk- 

k=l \fj.\=t-l~dcg{hk) k=l 0<|/j|<t-2-dog(/ifc) 

Since Xcis(2:f') is always a multiplicative variable for the polynomial {x^ /xcis(xt^))hk & 
T, we know that each polynomial in kerM(_i(y) can be represented by the poly- 
nomials in T and Ti, where 

Ti = {x,g \l<t< ch{g), geT, deg(.g) =t~2}. 

The polynomials in Ti and T have different leading terms, hence T U Ti is a 
linearly independent basis of ker Mt-i{y). Moreover, T is a reduced basis of 
kei Mt-2{y), and Ti consists of all linearly independent polynomials with degree 
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t — 1 in ker Mt~i{y). We can deduce that the number of polynomials in Ti is equal 
to corankMt_i(y) — corankMt_2(y). On the other hand, let aj denote the number 
of polynomials of class j and degree t — 2 in T. Since the set Ti is constructed 
by multiplying polynomials in T of degree t — 2 by their multiplicative variables 
only, the total number of polynomials in Ti is equal to JQ^i- Therefore, the 

condition ([TUl) is satisfied. □ 

3.4. An Extension to /(Vr(/) n A). Consider the semialgebraic set A := {x ^ 
\ fi{x)>0,..., fs{x) > 0} and the ^radical ideal /(Fr(/) HA). We restrict to 
a subset JCt^A ^ defined as 

ICt.A ICtn[ye M^- | Aft_d,„ (ry) h for all^/ e {0, 1}^} , 

where d/.^ = [deg(/'^)/2] and t>d = ma.xi<j<„i,i^e{os}={dj, d/-}. 

For t large enough, Lemma [276l and Theorem l2 . 7l show that the information about 
v^T is contained in the projection of a generic element y G ICt^A- Thus, propositions 
and theorems discussed above are true for generic elements y in JCt,A- 

The following theorem can be seen as a variant of Theorem 13.11 for the semi- 
algebraic set A. The proof uses exactly the same reason as in Theorem 13.91 and 
Theorem 13. 101 after replacing K-t and by K-t^A and \^ respectively. 

Theorem 3.11. Suppose the condition ilO\) holds for a generic element y G /Ct,yt, 
and t > 2d. Then a reduced basis of the null space of Mt-2{y) is a weak Pommaret 
basis of (ker Aft_2(y)) under the monomial ordering ^tdcg o-nd 

I C (kerMt_2(y)) C /(14(/) H A). 

Remark 3.12. For computing a Pommaret basis of (kerMt_2(y)), we add the 
defining polynomials {fi, ■ ■ ■ , fs} of the semialgebraic set A to the input of the 
above algorithm and additional constraints Mt-df^ (/"y) h for all i' e {0, 1}* to 
the semidefinite program (jl2p . 

4. Numerical examples 

We present here the results obtained by applying Algorithm l3.3l to some examples 
in [2711301 m [36] and others. 

Example 4.1. Consider the 2-dimensional ideal / = (/ii, /i2, /is) taken from |36l 
p.397, Eq. (9.60)] where 

hi = xl+ XiX2 - X1X3 - Xi- X2 + X3, 
/l2 = X1X2 + xl - X2X3 - Xi- X2+ X3, 
/l3 = X1X3 + X2X3 ~ xl - Xi - X2 + X3. 

The rank and corank sequences for truncated moment matrices Mt-i{y) are shown 
in Table 1 and 2. We set r — 10^^ and xi ^tdog X2 ^tdcg X3. For t=4, we have 

3 

jaj — 6, and corank M4_i — corank M4_2 = 6. 
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Hence, the condition ([TOl) is satisfied. The Pominaret basis computed by Algorithm 
[Qfor t = 4 is 

{Xi + X2 - - xl ~ X1X2 + XiX3,Xi + X2- X3 - X1X2 - xl+ X2X3, 

3xi + 3x2 — 3x3 ~ ^1 ~ 2xiX2 — X2 + X3}. 

From Table 3, we note that the condition ([TO)) is also satisfied for t = 5,6, 7. For this 
example, we can show that (ker M4-2 (y)) = V^, and a reduced basis of ker Af4_2 (y) 
is a Pommaret basis of -^7. Hence, the condition (fTO)) can be satisfied by arbitrary 
t > 4. 



Table 1. The rank of Mt-e{y) TABLE 2. The corank of Mt~i{y) 



Order 


e = o 


e = 1 


£^2 


Order 


£ = 


£=1 


£^2 


t=4 


16 


11 


7 


t=4 


19 


9 


3 


t=5 


22 


16 


11 


t=5 


34 


19 


9 


t=6 


29 


22 


16 


t=6 


55 


34 


19 


t=7 


37 


29 


22 


t=7 


83 


55 


34 



Table 3. The aj of a reduced basis of kerMt_2(y) 



Order ai a2 03 



t=4 


1 


1 


1 


t=5 


3 


2 


1 


t=6 


6 


3 


1 


t=7 


10 


4 


1 



1x1+1x2+1x3=6 
3x1+2x2+1x3=10 
6x1+3x2+1x3=15 
10x1+4x2+1x3=21 



Example 4.2. Consider the polynomial system P — {/ii,/i2} in ^30, p. 20, Ex 
1.4.6], where 

hi = xl- X2, 

h2 = X1X2 ~ X3. 

For the term order 0:3 ^tdog Xi ^tdcg 2^2, we have cls(a;i) = 2, cls(x2) = 3, cls(a;3) = 
1. Let T = 10~^, the rank and corank sequences for truncated moment matrices 
Mt-e{y) are shown in Table 4 and 5. 



Table 4. The rank of Mt-i{y) TABLE 5. The corank of Mt-i{y) 



Order 


^ = 


£ = l 


£ = 2 


Order 


£ = 


£ = 1 


^ = 2 


t=3 


12 


7 


4 


t=3 


8 


3 





t=4 


16 


10 


7 


t=4 


19 


10 


3 


t=5 


20 


13 


10 


t=5 


36 


22 


10 
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Table 6. The aj of a reduced basis of ker Aft_2(y) 



Order 




a2 


as 






t=3 











0x2+0x3+0x1= 


=0 


t=4 


2 


1 





2x2+1x3+0x1= 


=7 


t=5 


3 


1 


3 


3x2+1x3+3x1= 


12 



For t=4, we have 

3 

jaj = 7, and corankM4_i — corankAf4_2 = 7. 

Hence, the condition PU)) is satisfied. The Pommaret basis computed by Algorithm 
lOfor t = 4 is 

{xj - a;2, X1X2 - 2:3, xl - X1X3}. 
Example 4.3. Consider the ideal / = (/ii,/i2) in ^71 p.l23, Ex 7.41] with 

hi = xl + x\ + x\ — 2^ 
/12 — x\ + x\ — X:^. 

The real variety V]r(/) for this ideal is strictly contained in Vc{I)- We set r — 
10~^ and Xl ^tdcg X2 ^tdcg 2^3- The rank and corank sequences for truncated 
moment matrices Mt^iijj) are shown in Table 7 and 8. 



Table 7. The rank of Mt-e{y) TABLE 8. The corank of Mt-i{y) 



Order 




e = i 


£ = 2 


Order 


^ = 


£=1 


^ = 2 


t=3 


7 


5 


3 


t=3 


13 


5 


1 


t=4 


9 


7 


5 


t=4 


26 


13 


5 


t=5 


11 


9 


7 


t=5 


45 


26 


13 


t=6 


13 


11 


9 


t=6 


71 


45 


26 



Table 9. The aj of a reduced basis of kerMt_2(y) 



Order 


ai 


a2 


as 










t=3 








1 


0x1- 


f0x2 


4-1x3= 


=3 


t=4 


1 


2 


1 


1x1- 


f2x2- 


4-1x3= 


=8 


t=5 


4 


3 


1 


4xH 


-3x2H 


-1x3= 


+3 


t=6 


8 


4 


1 


8x14 


-4x2H 


-1x3= 


+9 



For t=4, we have 
3 

jckj = 8, and corank M4_i — corank M4_2 = 8. 
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Hence, the condition ([TOl) is satisfied. The Pommaret basis computed by Algorithm 
[Qfor t = 4 is 

For this example, we can show that (kerM4_2(y)) — v^, and a reduced basis of 
kerM4_2(y) is a Pommaret basis of v^- Hence, the condition ([TU]) can be satisfied 
by arbitrary t > 4. 

Example 4.4. Consider the ideal / — (/ii, /12, /is) in [351 p. 61, Ex 2.4.12], where 

hi = xl+ X2X3 - x\, 
h2 = xiX'i + a;ia;2 - X3, 
/13 = X2XZ + xl + xl — Xi. 

Let T — lO"*" and the term order be xi ^tdcg X2 ^tdcg X3. The rank and corank 
sequences for the truncated moment matrices Mt-i{y) are shown in Table 10 and 
11. 



Table 10. The rank of Mt-t{y) TABLE 11. The corank of Mt-e{y) 



Order 




e = i 


e = 2 


Order 




£ = 1 


i = 2 


t=4 


13 


10 


7 


t=4 


22 


10 


3 


t=5 


16 


13 


10 


t=5 


40 


22 


10 


t=6 


19 


16 


13 


t=6 


65 


40 


22 


t=7 


22 


19 


16 


t=7 


98 


65 


40 



Table 12. The aj of a reduced basis of kerMt-2(y) 



Order 




a2 


as 




Ej=i jo^j 




t=4 


1 


1 


1 


1x1- 


flx2+lx3= 


=6 


t=5 


4 


2 


1 


4xH 


-2x2+1x3= 


11 


t=6 


8 


3 


1 


8xU 


-3x2+1x3= 


17 


t=7 


13 


4 


1 


13x1- 


f4x2+lx3= 


=24 



The condition ([TOl) can not be satisfied for t from 4 to 7. Actually, Seller showed 
in |32j that the coordinates {xi,X2,X3) are not (5-regular for the ideal /. However, 
if we perform the linear transformation suggested in j32| . xi = X3, £2 = X2 + x^, 
X3 — xi, after an auto-reduction, we obtain the polynomial system P = {xiX2 — 
i§, ai2i3 — ii, ^2 ~ xs}- We choose an ordering xi ^tdog X2 ^tdcg is and r = 10^*. 
The rank and corank sequences for the truncated moment matrices Mt-t{y) are 
shown in Table 13 and 14. 

For t=4, we have 

3 

jaj — 7, and corank M4_i — corank M4_2 = 7. 
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Table 13. The rank of Mt^e{y) TABLE 14. The corank of Mt-e(y) 



Order 


£ = £ = 1 


^ = 2 


Order £ = l^l 


£ = 2 


t=4 


13 10 


7 


t=4 22 10 


3 


t=5 


16 13 


10 


t=5 40 22 


10 


t=6 


19 16 


13 


t=6 65 40 


22 




Table 15. 


The a-j of 


a reduced basis of kerMt-2(y) 






Order 


cti a2 








t=4 


2 


1 0x1+2x2+1x3=7 






t=5 


3 3 


1 3x1+3x2+1x3=12 






t=6 


7 4 


1 7x1+4x2+1x3=18 





Hence, the condition (flUl) is satisfied. The Ponimaret basis computed by Algorithm 
[Qfor t = 4 is 

{-a;3 + xl,-xi + X2X'i, -X1X2 + xj}. 

Example 4.5. Consider the ideal / = {hi, /12), where 

hi = {xi ~ X2)ixi + X2)'^{xi + x\+ 2:2), 

h2 = (a;i - a;2)(a;i + X2Y{xl + x^). 

In this example, / is not a radical ideal. We set r = 10^** and xi ^tdcg X2- The 
rank and corank sequences for the truncated moment matrices Mt-i{y) are shown 
in Table 16 and 17. 

Table 16. The rank of Mt-t(y) TABLE 17. The corank of Mt-t(y) 



Order 




t = 1 


£ = 2 


Order 


^ = 


£= 1 


£ = 2 


t=7 


15 


13 


11 


t=7 


21 


15 


10 


t=8 


17 


15 


13 


t=8 


28 


21 


15 


t=9 


19 


17 


15 


t=9 


36 


28 


21 



Table 18. The qj of a reduced basis of kerMt-2(t/) 



Order 




a2 




=1 


t=7 


3 


1 


3x1- 


flx2=5 


t=8 


4 


1 


4x1- 


flx2=6 


t=9 


5 


1 


5x1- 


flx2=7 



For t=7, we have 
2 

j'ttj = 5, and corank My-i — corank M7_2 ~ 5. 
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Hence, the condition ([TOl) is satisfied. The Pominaret basis computed by Algorithm 
[Qfor t = 7 is 

{-xl + xl}. 

It should be noticed that for this example, if we set tolerance t < 10^^, the rank and 
corank sequences for the truncated moment matrices Mt-eiu) will be completely 
different from those shown in Table 16 and 17, and we can not get {—xf + x\} as 
a Pommaret basis of v^. 

Example 4.6. We compute I{Vs,[I) n A) for / ^ (/ii, /i2), 

hi = {xi ~ X2){xi + X2){xi +x\+ X2), 
h2 = (xi - X2)ixi + X2){xl + xl), 

and 

A = {{XUX2) e I > 1, a;2 > 1}. 

Let us set r — 10^^ and xi ^tdcg ^2, the rank and corank sequences for the 
truncated moment matrices Mt-i(ij) with y g ICt.A shown in Table 19 and 20. 



Table 19. The rank of Mt-eiy) TABLE 20. The corank of Mt-e{y) 



Order 


£ = e=i 


1^2 


Order 


£ = 


£^l 


£^2 


t=6 


8 6 


5 


t=6 


20 


15 


10 


t=7 


9 7 


6 


t==7 


27 


21 


15 


1=8 


10 8 


7 


t=8 


35 


28 


21 



Table 21. The qj of a reduced basis of ker Mt-2{y) 



Order 


ai 


a2 




Li jaj 




t=6 


3 


1 


3x1- 


flx2= 


:5 


t=7 


4 


1 


4x1- 


f 1x2= 


6 


t=8 


5 


1 


5x1- 


flx2= 


-7 



For t=6, we have 
2 

^^j'ttj = 5, and corank A^e_i — corank M6_2 = 5. 

Hence, the condition (|10p is satisfied. The Pommaret basis we obtain by Algorithm 
lOfor i = 6 is 

{-Xl + X2} 

for I{VRiI)nA). 



a certificate for semidefinite relaxations in computing positive dimensional real varieties 

5. Conclusion 

In this paper we present a semidefinite characterization for computing a Pom- 
maret basis of an ideal J, where J is generated by polynomials in the kernel of a 
truncated moment matrix and satisfies I Q J Q /(Vk(7)). Our approach is stim- 
ulated by the previous work inJTIlIIHlESllMlEIlEIllMlISS]. By combining the 
geometric involutive theory with the results on positive semidefinite moment matri- 
ces, we introduce a new stopping condition (|10|) for the semidefinite program p2|) 
and prove the finite termination of the algorithm in a (5-regular coordinate system. 
Although from the tables in Section |4l we can check that the condition ((T0|) can be 
satisfied by higher order moment matrices once it is satisfied at some order, in gen- 
eral, we can not guarantee this property. Therefore, unlike flat extension conditions 
proposed by Curto and Fialkow in [8] for finite rank moment matrices, we can not 
prove the computed Pommaret basis is an involutive basis of the real radical ideal. 
Finally, we wish to mention that results computed by semidefinite programming 
and numerical linear algebra are approximate. Therefore, our condition (jlOl) can 
only be checked with respect to a given tolerance. For improperly chosen tolerance, 
we might not be able to give a meaningful answer. 
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